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ABSTRACT 

We study numerically the dependence of the critical magnetic Reynolds number Rm c for the turbulent small- 
scale dynamo on the hydrodynamic Reynolds number Re. The turbulence is statistically homogeneous, isotropic, 
and mirror-symmetric. We are interested in the regime of low magnetic Prandtl number Pm = Rm/Re < 1, 
which is relevant for stellar convective zones, protostellar disks, and laboratory liquid-metal experiments. The 
two asymptotic possibilities are Rm c — * const as Re — ► oo (a small-scale dynamo exists at low Pm) or Rm c /Re = 
Pm c — > const as Re — > oo (no small-scale dynamo exists at low Pm). Results obtained in two independent sets of 
simulations of MHD turbulence using grid and spectral codes are brought together and found to be in quantitative 
agreement. We find that at currently accessible resolutions, Rm c grows with Re with no sign of approaching a 
constant limit. We reach the maximum values of Rm c ~ 500 for Re ~ 3000. By comparing simulations with 
Laplacian viscosity, fourth-, sixth-, and eighth-order hyperviscosity and Smagorinsky large-eddy viscosity, we 
find that Rm c is not sensitive to the particular form of the viscous cutoff. This work represents a significant 
extension of the studies previously published by Schekochihin et al. (2004a) and Haugen et al. (2004a) and the 
first detailed scan of the numerically accessible part of the stability curve Rm c (Re). 

Subject headings: magnetic fields — methods: numerical — MHD — turbulence 



The magnetic Prandtl number Pm, which is the ratio of the 
kinematic viscosity to the magnetic diffusivity, is a key pa- 
rameter of MHD turbulence. In fully ionized plasmas, Pm w 
2.6 x 10~ 5 T 4 /n, where T is the temperature in kelvins and 
n the ion number density in units of cm -3 . In hot rarefied 
plasmas, such as the warm and hot phases of the interstel- 
lar medium or the intracluster medium, Pm 3> 1. In contrast, 
in the Sun's convective zone, Pm ~ 10~ 7 to 10 -4 , in planets, 
Pm ~ 10~ 5 , and in protostellar disks, while estimates vary, it is 
also believed that Pm <C 1 (e.g., Brandenburg & Subramanian 
2004). All these astrophysical bodies have disordered fluctuat- 
ing small-scale magnetic fields and, in some cases, also large- 
scale "mean" fields. As they also have large Reynolds numbers 
and large-scale sources of energy, they are expected to be in a 
turbulent state. It is then natural to ask if their magnetic fields 
are a product of turbulent dynamo. 

To be precise, there are two types of dynamo. The large- 
scale, or mean-field dynamo generates magnetic fields at scales 
larger than the energy-containing scale of the turbulence, as is, 
for example, the case in helical turbulence. The small-scale dy- 
namo amplifies magnetic fluctuation energy below the energy- 
containing scale of the turbulence. The small-scale dynamo is 
due to random stretching of the magnetic field by turbulent mo- 
tions and does not depend on the presence of helicity. Mean- 
field dynamos typically predict field growth on time scales as- 
sociated with the energy-containing scale (or longer), while the 
small-scale dynamo amplifies magnetic energy at the turbulent 
rate of stretching. Thus, the small-scale dynamo is usually a 
much faster process than the mean-field dynamo, and the large- 



scale field produced by the latter can be treated as approxi- 
mately constant on the time scale of the small-scale dynamo. 
The mean-field dynamo (or, more generally, a large-scale mag- 
netic field of any origin) also gives rise to small-scale magnetic 
fluctuations because of the turbulent shredding of the mean 
field: this leads to algebraic-in-time growth of the small-scale 
magnetic energy — again, a slower generation process than the 
exponential-in-time small-scale dynamo. 

In the systems with Pm 3> 1, the existence of the small-scale 
dynamo is well established numerically and has a solid theoret- 
ical basis (see Schekochihin et al. 2004b for an account of the 
relevant theoretical and numerical results and for a long list of 
references). The situation is much less well understood for the 
case of small Pm. It has been largely assumed that a small-scale 
dynamo should be operative in this regime as well. For exam- 
ple, the presence of large amounts of small-scale magnetic flux 
in the solar photosphere (e.g., Title 2000) has been attributed 
to small-scale dynamo action. This appeared to be confirmed 
by numerical simulations of the MHD turbulence in the con- 
vective zone (Cattaneo 1999; Cattaneo et al. 2003; Nordlund 
2003). However, such simulations are usually done at Pm > 1 
(Pm = 5 in Cattaneo's simulations). Previous attempts to sim- 
ulate MHD turbulence in various contexts with Pm < 1 found 
achieving dynamo in this regime much more difficult than for 
Pm > 1 (Nordlund et al. 1992; Brandenburg et al. 1996; Nore et 
al. 1997; Christensen et al. 1999; Maron et al. 2004). A system- 
atic numerical investigation of the effect of Pm on the efficiency 
of the small-scale dynamo was carried out by Schekochihin et 
al. (2004a), who found that the critical magnetic Reynolds num- 
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ber Rm c required for the small-scale dynamo to work increases 
sharply at Pm < 1 . An independent numerical study by Haugen 
et al. (2004a) confirmed this result. 

What are the basic physical considerations that should guide 
us in interpreting this result? First of all, let us stress that all 
working numerical small-scale dynamos are of the large-Pm 
kind (the case of Pm = 1 is nonasymptotic, but its properties that 
emerge in numerical simulations suggest that it belongs to the 
same class). Two essential features of the large-Pm dynamos 
are (1) the scale of the velocity field is larger than the scale of 
the magnetic field, and (2) the velocity field that drives the dy- 
namo is spatially smooth and locally looks like a random linear 
shear, so the dynamo is due to exponential-in-time separation of 
Lagrangian trajectories and the consequent exponential stretch- 
ing of the magnetic field. The basic physical picture of such 
dynamos (Zeldovich et al. 1984; see discussion in Schekochi- 
hin et al. 2004b; see also a review of an alternative but com- 
plementary approach by Ott 1998) explicitly requires these two 
conditions to hold. The map dynamos and the dynamos in de- 
terministic chaotic flows that were extensively studied in the 
1980s-1990s (see review by Childress & Gilbert 1995) are all 
of this kind. The numerical dynamos with Pm > 1 (the first 
due to Meneguzzi et al. 1981) are of this kind as well because 
they are driven by the spatially smooth viscous-scale turbulent 
eddies, which have the largest turnover rate. 
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FIG. 1. — Sketch of the two possible shapes of the stability curve Rm c vs. Re 
for the small-scale dynamo. 

When Pm <C 1 with both Rm 3> 1 and Re ^> 1 , the character- 
istic scale Ib of the magnetic field lies in the inertial range. For 
Kolmogorov turbulence, a simple estimate gives l B ~ Rm~ 3 / 4 /o, 
where Iq is the energy-containing scale. As the viscous scale 
is l v ~ Re~ 3//4 /o, we have Iq <C Ib <C l v - In a rough way, one 
can think of the turbulent eddies at scales / > Ib as stretching 
the field at the rate u//l and of the eddies at scales I < l B as 
diffusing the field with the turbulent diffusivity m;Z. In Kol- 
mogorov turbulence, m/ ~ Z 1 / 3 , so both the dominant stretching 
and the dominant diffusion are due to the eddies at scale I ~ Is- 
The resulting rates of stretching and of turbulent diffusion are 



of the same order, so the outcome of their competition cannot 
be determined on this qualitative level (Rraichnan & Nagara- 
jan 1967). An important conclusion, however, can be drawn. 
If the bulk of the magnetic energy is at the scale l B , the ex- 
istence of the dynamo is entirely decided by the action of the 
velocities at the scale Is- Then it cannot matter where in the 
inertial range Ib lies. But ~ Pm~ 3 ^ 4 , so the value of Pm 
does not matter as long as it is asymptotically small. Therefore, 
there are two possibilities: either there is a dynamo at low Pm 
and Rm c — > const as Re — > oo or there is not and there exists 
a finite Pm c = Rm c /Re — > const as Re — > oo. Strictly speak- 
ing, the third possibility is that Rm c oc Re Q , where a is some 
fractional power, but this can only happen if the intermittency 
of the velocity field (non-self-similarity of the inertial range) is 
important for the existence of the dynamo. 9 

The two possibilities identified above are illustrated in Fig. 1 . 
Several authors (Vainshtein 1982; Rogachevskii & Kleeorin 
1997; Boldyrev & Cattaneo 2004) showed that, given certain 
reasonable assumptions, the first possibility (Rm c — > const) is 
favored by the Kazantsev (1968) model: the small-scale dy- 
namo in a Gaussian white-in-time velocity field. In particular, 
Boldyrev & Cattaneo (2004) found that the Kazantsev model 
gives Rm c that is roughly 10 times larger in the Pm <C 1 regime 
than in the Pm 3> 1 regime (Rogachevskii & Kleeorin 1997 pre- 
dict Rm c ~ 400, which is consistent with that). This prompted 
them to declare the issue settled on the grounds that the failure 
of the dynamo in numerical experiments at current limited res- 
olutions is compatible with such an increase in Rm c . However, 
the Pm <C 1 dynamo in the Kazantsev model is a quantitative 
mathematical consequence of the model, and it is not known if 
and how it is affected by such drastic and certainly unrealistic 
assumptions as the Gaussian white-noise statistics for the ve- 
locity field. 10 The existence of a dynamo in real turbulence is 
also a quantitative question (see discussion above), so it cannot 
be decided by a model that is not a quantitative approximation 
of turbulence. 

Thus, the issue cannot be considered settled until definitive 
numerical evidence is produced. This is an especially hard task 
because we do not know just how high a magnetic Reynolds 
number we must achieve in order to clearly see the distinction 
between Rm c — > const and Rm c /Re — » const. In this Letter, we 
have collected numerical results from two independent compu- 
tational efforts: simulations using an incompressible spectral 
MHD code (see code description in Maron & Goldreich 2001; 
Maron et al. 2004) and weakly compressible simulations using 
a grid-based high-order MHD code (the PENCIL CODE 11 ). 

The equations we solved numerically (in a triply periodic 
cube) are 

a,u + u.Vu = -^ + ^4^ + F vlsc+ f, (1) 
p 4-irp 

9,B = Vx(uxB) + i)V 2 B, (2) 



where u is the velocity and B is the magnetic field (the PEN- 
CIL CODE, in fact, solves the evolution equation for the vector 
potential A and then computes B = V x A). All runs reported 
below are in the kinematic regime, |B| <C |u|, so the Lorentz 
force in Eq. (1) plays no role. Turbulence is driven by a random 

9 The role of coherent structures can be prominent in quasi-two-dimensional dynamos (three-component velocity field depending on two spatial variables), where the 
inverse cascade characteristic of the two-dimensional turbulence gives rise to persistent large-scale vortices, which drive the dynamo (Smith & Tobias 2004). 

10 Vainshtein & Kichatinov (1986) argue that the equations that arise from the Kazantsev model are valid for non-white velocity fields if n-point joint probability 
density functions of Lagrangian displacements satisfy Fokker-Planck equations with some diffusion tensor. They further assume (on dimensional grounds) that 
this diffusion tensor scales as the scale-dependent turbulent diffusion ~ u\L This is, in fact, a closure scheme that we believe to be equivalent to Kazantsev's 
zero-correlation-time theory. 

1 1 See http://www.nordita.dk/software/pencil-code. 
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white-in-time nonhelical body force f concentrated at k = ko, 
where ko is the wavenumber associated with the box size. The 
(hyper)viscous force is 



F vhc =Iv-[2(p)i/„(-V 2 )"" 1 S] 1 
p 



where v n is the fluid viscosity and 



1 / dui duj\ 1 
^ = 2^ + &-J-3^ V ' U - 



(3) 



(4) 



In the spectral simulations, the density p = 1, and the incom- 
pressibility constraint V • u = is enforced exactly via the deter- 
mination of the pressure p. The grid simulations are isothermal: 
p = c\p with sound speed c s = 1, and the density satisfies 



a,p+v-(pu) = o. 



(5) 



We stay in the weakly compressible regime of low Mach num- 

l li 

bers M = (u 2 ) /c s ~ 10 -1 and p ~ (p) = 1 (angular brackets 
denote volume averages). Some numerical results on the onset 
of dynamo action at larger Mach numbers are given in Haugen 
et al. (2004b). 

The dissipation in the induction equation (2) is always Lapla- 
cian with magnetic diffusivity T) (we choose not to tamper with 
magnetic dissipation because we are interested in the sensitive 
question of field growth or decay). With regard to the viscous 
dissipation, we perform three kinds of simulations: 

1. Simulations with Laplacian viscosity: n = 1 in Eq. (3). 

2. Simulations with fourth-, sixth-, and eighth-order hyper- 
viscosities: n = 2, 3, and 4, respectively, in Eq. (3). 

3. Large-eddy simulations (LES) with the Smagorinsky ef- 
fective viscosity (e.g., Pope 2000): in Eq. (3), n = 1, and v\ is 

replaced by v$ = (Cs A) (2S : S) , where A is the mesh size 
and Cs = 0.2 is an empirical coefficient. 

The magnetic Reynolds number is defined by Rm = 

1 /2 

(u 2 ) /korj, where ko is the box wavenumber (the smallest 
wavenumber in the problem). For the runs with Laplacian vis- 
cosity (n = 1), the hydrodynamic Reynolds number is Re = 

1/2 

(u 2 ) /hjV\, For hyperviscous runs and for LES, we define 
Re by replacing v\ with the effective viscosity: 



^eff = e/(2S:S)=e/(|Vu| 2 } 



(6) 



(the second expression is for the spectral simulations, where 
V • u = exactly). Here e = (f ■ u) is the total injected power 
and is equal to the total energy dissipation. As the forcing f 
is white in time, e = const: indeed, given (/'(? ,x)fj(t' ,x')) = 
(5(f-f')e' 7 (x-x'), it is easy to show that e = \ e"(0). 

The results of all our simulations are presented in Fig. 2, 
where Rm c is plotted versus Re. Each value of Rm c 
was computed by interpolating between least-squares-fitted 
growth/decay rates of a growing and a decaying run. Error bars 
are based on Rm and Re for these pairs of runs. The only ex- 
ception is the point enclosed in a circle, which corresponds to 
(Rm, Re) for a run that appeared to be marginal (in this case we 
could not afford the resolution necessary to achieve a decay- 
ing case). The run times in all cases were long enough for the 
least-squares-fitted growth rates to stop changing appreciably 
(typically this required about 20 box-crossing times, but cases 
close to marginal needed longer running times). 



• Laplacian (JLM) 

□ <1t.h order hyper ('JT ,M) 
6th — order hyper/(JLM) 

- O 8th-order hyper (JLM) 

A Laplacian (PENCIL) 

-X- 6th -order hyper (PENCIL) 

h A Smagorinsky (PENCIL) 





Pm=10 



Pm=l 



Pm=0.1 



Pm = 0.01 



10 z 
Rc 



10 3 



10* 



FIG. 2. — Dependence of Rm c on Re. "JLM" refers to simulations done 
with the incompressible spectral code written by J. L. Maron: runs with Lapla- 
cian viscosity, fourth-, sixth-, and eighth-order hyperviscosity (resolutions 64 3 - 
256 3 ). In this set of simulations, hyperviscous runs were done at the same val- 
ues of i) as the Laplacian runs, so the difference between the results for these 
runs is nearly imperceptible. "PENCIL" refers to weakly compressible simula- 
tions done with the PENCIL CODE: runs with Laplacian viscosity, sixth-order 
hyperviscosity, and Smagorinsky large-eddy viscosity (resolutions 64 3 -512 3 ). 

We see that there is good agreement between the results for 
runs with different forms of viscous dissipation; this confirms 
the natural assumption that the field-generation properties of 
the turbulence at low Pm are not sensitive to the way the veloc- 
ity spectrum is cut off. It is also encouraging that results from 
two very different codes are in quantitative agreement. 

Our previous studies (Schekochihin et al. 2004a; Haugen et 
al. 2004a) had the maximum value of Rm c ~ 200. The results 
reported here raise it to ~ 500, with the corresponding values of 
Pm c around 0.15. While a roughly 10-fold increase with respect 
to Rm c for the Pm = 1 dynamo has now been achieved, there is 
thus far no sign of Rm c reaching an asymptotically constant 
value. This said, the current resolutions are clearly still insuffi- 
cient to make a definitive judgement, although we are now very 
close to values of Rm c predicted by the theories based on the 
Kazantsev model (Rogachevskii & Kleeorin 1997; Boldyrev & 
Cattaneo 2004) — whether or not the model yields quantita- 
tively correct predictions should become clear in the near fu- 
ture. 

The numerical results reported above concerned the depen- 
dence Rm c (Re) for the turbulent small-scale dynamo, i.e., the 
ability of turbulent velocity fluctuations to amplify magnetic 
energy at scales smaller than the energy-containing scale of the 
turbulence. The Rm c (Re) dependence is also an interesting is- 
sue for other kinds of dynamo. 

If the velocity field is non-mirror-symmetric, it can often 
drive the mean-field dynamo (MFD), which means the growth 
of the magnetic field at scales larger than the energy-containing 
scale of the turbulence (Krause & Radler 1980). This large- 
scale field generated by the MFD, just like a mean field imposed 
externally, can induce small-scale magnetic fluctuations as it is 
shredded by the turbulence, so the total field has both a mean 



4 



SCHEKOCHIHIN ET AL. 



(large-scale) and a fluctuating component. In many cases, the 
breaking of the mirror symmetry leads to a non-zero value of 
the net helicity (u • ( V x u)) ^ (the average is over all scales 
smaller than the energy-containing scale of the turbulence, not 
over the entire volume of the system). The mean-field gen- 
eration is then referred to as the a-effect. The stability curve 
Rm c (Re) for the a-effect is different than for the small-scale dy- 
namo: it is essentially a condition for at least one unstable large- 
scale mode to fit into the system. In a numerical study done 
with the same code as the grid simulations reported above but 
with fully helical random forcing, Brandenburg (2001) found 
much lower values of Rm c than for the small-scale dynamo and 
very little dependence of Rm c on Pm for Pm > 0.1. 

A non-zero net helicity is not a necessary condition for the 
MFD (e.g., Gilbert et al. 1988). In fact, it has been suggested 
recently by Rogachevskii & Kleeorin (2003) that the MFD can 
be driven simply by the presence of a constant mean velocity 
shear (shear-current or (5-effect) — a very generic possibility 
of obvious relevance to systems with mean flows. Mean flows 
are present in many astrophysical cases and in all current lab- 
oratory dynamo experiments (Gailitis et al. 2004; Miiller et al. 
2004; Bourgoin et al. 2002; Lathrop et al. 2001; Forest et al. 
2002). A mean flow can be a dynamo in its own right: an MFD 
(field growth at scales above the flow scale) and, if the flow has 
chaotic trajectories in three dimensions, also a small-scale dy- 
namo (field growth at scales ~ Rm~ 1,/2 times the scale of the 
flow; see Childress & Gilbert 1995 — as noted above, small- 
scale dynamos in deterministic chaotic flows are equivalent to 
the large-Pm case). When Re is large, the energy of the tur- 
bulent velocity fluctuations is comparable to the energy of the 
mean flow. The critical Rm required for field growth will have 
some dependence on Re, which reflects the effect of the turbu- 
lence on the structure of the mean flow and/or on the effective 
value of the magnetic diffusivity (the /3-effect; see Krause & 
Radler 1980). This dependence was the subject of two recent 



numerical studies: of the dynamo in a turbulence with a con- 
stant Taylor-Green forcing by Ponty et al. (2005), and of the 
Madison dynamo experiment (propeller driving in a spherical 
domain) by Bayliss & Forest (2004). The Re dependence of 
Rm c that emerges from such simulations is distinct from that 
for a pure small-scale dynamo. Indeed, Y. Ponty et al. (2005, 
private communication) have shown that, in the limit of large 
Re, the value of Rm c in their simulations tends to a constant 
that coincides with Rm c calculated for the mean flow alone, i.e., 
for the velocity field with fluctuations removed by time averag- 
ing. In contrast, the subject of the present Letter has been the 
possibility of a small-scale dynamo driven solely by turbulent 
fluctuations, in the absence of a mean flow. The importance of 
this possibility or lack thereof is that such a dynamo, if it ex- 
ists, occurs at the turbulent stretching rate associated with the 
resistive scale. This is much faster (by a factor of ~ Rm'^ 2 ; 
see discussion above) than the growth rate of any MFD or of 
a small-scale dynamo associated with the mean flow, i.e., than 
the stretching rate at the energy-containing scale or at the scale 
of the mean flow. 
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